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Calculations  of  Dielectric  Constant  for  AlGalnAs  Quaternary  Semiconductor  Alloy  in  the 
Transparent  Region  and  Above  (0.4-4.0eV) 

M.  Linnik  and  A.  Christou 

Department  of  Materials  and  Nuclear  Engineering  and  Materials  Research  Science  and 
Engineering  Center,  University  of  Maryland,  College  Park  MD  20742 

ABSTRACT 

The  modeling  of  the  spectral  behavior  of  the  refractive  index  of  AlGalnAs  quaternary  III- 
V semiconductor  alloy  in  the  energy  range  from  0.4  to  4eV,  including  the  transparent  region,  is 
presented.  The  extended  model  of  interband  transition  contributions  incorporates  not  only  the 
fundamental  absorption  edge  contribution  to  the  dielectric  function,  but  also  contributions  from 
higher  energy  and  indirect  transitions.  It  is  demonstrated  that  indirect  energy  transitions  must  be 
included  in  the  calculations  of  the  complex  dielectric  function  of  the  material  in  the  transparent 
region.  Indirect  transitions  from  different  critical  points  in  the  Brillouin  zone  are  treated 
separately.  The  comparison  between  the  theoretical  refractive  indices  and  the  experimental  data 
for  AlGalnAs  alloy  is  presented.  These  calculations  have  been  applied  to  the  design  of  Bragg 
mirrors  with  the  highest  refractive  index  contrast  for  heterostructure  lasers. 

INTRODUCTION 

The  design  and  analysis  of  such  devices  as  injection  lasers,  photodiodes,  detectors,  solar 
cells,  multilayer  structures,  and  microcavities  requires  the  exact  knowledge  of  the  optical 
constants  of  III-V  compound  semiconductors  in  the  region  near  the  fundamental  absorption  edge 
as  well  as  at  the  higher  photon  energies.  In  modeling  of  the  optical  constants  of  semiconductors 
in  the  fundamental  optical  region,  several  approaches  are  typically  used:  (1)  empirical  formulas, 
(2)  damped  harmonic  oscillator  (DHO)  models,  (3)  standard  critical  point  (SCP)  models.  Optical 
constants  determined  from  empirical  formulas  (such  as  the  Sellmeier  dispersion  equations  for  the 
refractive  index  and  Urbach’s  rule  for  the  absorption  coefficient  [1],  or  the  expression  for  n 
based  on  interpolation  of  a dielectric  quantity  using  Vegard’s  rule  by  Burkhard  et  al.[ 2])  are  not 
related  through  the  Kramers-Kronig  dispersion  relation  and  are  valid  only  over  a very  limited 
energy  range. 

A semi-empirical  single  effective  oscillator  model  based  on  quasi-classical  Boltzmann 
equation  or  Drude  theory  proposed  by  Wemple  et  al.  [3]  does  provide  an  analytical  expression  for 
the  dispersion  of  the  semiconductor  refractive  index  at  photon  energies  significantly  below  the 
direct  band  edge.  The  Drude  theory  ignores  the  carrier  related  effects  around  the  band  gap,  and 
thus  the  results  are  valid  only  in  the  low  optical  frequency  region.  This  model  also  lacks  the 
agreement  with  experimental  data  at  the  band  edge,  which  is  the  energy  range  of  the  most 
interest  for  semiconductor  laser  devices. 

The  standard  critical  point  (SCP)  model  can  determine  the  position  of  critical  points  of 
the  semiconductor  band  structure,  but  cannot  accurately  predict  the  dielectric  function  [4],  The 
modified  SCP  model  was  initially  proposed  by  Korovin  [5]  and  Cardona  et  al.[ 6],  and  then 
developed  by  Adachi  [7],  and  Lin  et  cz/.[8].  The  model  of  interband  transition  contributions  (ITC 
model)  was  introduced  as  a method  to  analyze  the  refractive  index  of  III-V  compounds  at 
energies  below  and  above  the  direct  band  gap  by  including  the  electron-hole  pair  transitions,  and 
by  adding  the  excitonic  terms  at  the  two  lowest  energy  gap  transitions.  The  comparison  between 
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available  experimental  results  of  the  spectral  behavior  of  III-V  compound  semiconductors  and 
the  theoretical  data  calculated  using  the  above  mentioned  models  often  reveals  a lack  of 
agreement,  which  is  pronounced  for  the  photon  energies  around  the  fundamental  absorption 
edge.  These  differences  may  arise  from  the  excitonic  effects,  which  are  largely  ignored  in  the 
calculations  of  the  real  part  of  the  dielectric  constant  [8], 

In  the  present  work,  an  extended  model  of  interband  transition  contributions  (EITC)  is 
developed  for  the  calculations  of  real  and  imaginary  parts  of  the  dielectric  constant  of  compound 
semiconductors.  The  model  introduces  (1)  the  broadening  effects,  caused  by  phonon  and  defect 
scattering  in  direct  and  indirect  transitions;  (2)  the  strength  of  direct  band  gap  transitions  as  a 
function  of  the  effective  electron,  heavy  hole  and  light  hole  masses  of  the  semiconductor;  (3)  the 
exciton  contributions;  (4)  the  separate  contributions  of  Ex  and  Ei  indirect  band  gap  transitions  to 
the  real  and  imaginary  part  of  the  dielectric  constant.  The  importance  of  indirect  and  higher 
direct  energy  transitions  is  demonstrated  through  these  calculations  and  through  the  comparison 
with  experimental  results.  The  detailed  description  of  our  extended  ITC  model  is  given  in 
Reference  [9],  where  the  index  of  refraction  for  the  following  alloys  has  also  been  calculated  and 
compared  with  the  available  experimental  results:  A1P,  AlAs,  AlSb,  GaP,  GaAs,  GaSb,  InP, 

InAs,  InSb,  AIGaAsSb,  AlGalnAs,  AlGalnP,  GalnAsSb,  and  GalnPAs.  In  this  work  we  apply 
the  extended  ITC  model  for  the  calculation  of  the  dielectric  constant  of  AlGalnAs  alloy,  which  is 
then  used  for  the  Bragg  mirror  design. 

THE  EXTENDED  ITC  MODEL 

The  dielectric  constant  e(E)  = edEJ+fclE)  describes  the  optical  response  of  the  medium 
as  a function  of  photon  energy  E.  The  imaginary  part  of  the  dielectric  function  62(E)  is  calculated 
based  on  a simplified  model  of  the  band  structure  using  the  joint  density  of  states  for  each 
Critical  Point  (CP)  considered.  The  real  part  of  the  dielectric  function  6i(E)  was  calculated 
through  the  knowledge  of  the  imaginary  part,  62(E),  by  employing  the  Kramcrs-Kronig  relation 

[10] .  Thus,  the  total  imaginary  and  real  parts  of  the  dielectric  function  are  presented  as  a sum  of 
several  terms  that  represent  the  contribution  of  different  energy  CPs.  These  points  are  associated 
with  electronic  transitions  in  the  band  structure  at  the  energies  designated  as  E(l,  Eo+A,  En'\  Ei, 
Ej,  and  Ej.  In  case  of  quaternary,  AxByCzD,  semiconductor  alloy  each  of  the  terms  become  a 
function  of  the  alloy  mole  fraction,  x,  y and  z=l  -x-y.  There  arc  several  absorption  mechanisms 

[11]  that  contribute  to  the  imaginary  part  of  the  dielectric  constant,  therefore  62(E)  can  be  written 
as: 

£,(  E.x.y  ) = £■?'(  E,x,y ) + £"r"(  E,x,y ) + £*"'  J(  E,x, y ) + 

+ £_f'(  E.x.y ) + £?2(  E.x.y ) + £?(  E.x.y ) 

where  62^  and  62Eo+A  are  contributions  due  to  the  absorption  by  direct  interband  optical 
transitions  near  the  fundamental  absorption  edge  and  spin-orbit  transitions,  62™°  is  due  to  the 
absorption  by  the  discrete  series  of  excitons  near  the  En  energy  gap,  62EI  and  £2E2  are 
contributions  of  the  higher  energy  interband  transitions,  and  62 El  is  due  to  the  indirect  interband 
absorption  effects.  Similarly,  the  real  part  of  the  dielectric  function  can  be  presented  as  the 
following  sum: 

£,(  E.x.y ) = ef'(  E.x.y ) + ef"(  E.x.y ) + ef"' '(E.x.y ) + 

+ £,'(  E.x.y  ) + £,'-(  E.x.y  ) + ef‘(  E.x.y ) 
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The  absorption  spectrum  for  photon  energies  greater  than  the  band  gap  energy  is 
composed  of  many  peaks  correlated  with  Van-Hove  singularities  of  the  joint  density  of  states 
[1 1],  For  III-V  zinc-blende  type  semiconductors,  the  contributions  of  two  main  peaks  (Ei  and  E2) 
must  also  be  included.  These  peaks  correspond  to  the  direct  optical  transitions  at  the  L and  X 
points  of  the  BZ,  respectively.  The  E]  peak  is  treated  as  a two-dimensional  Mo  type  critical  point, 
while  the  structure  of  the  E2  peak  is  characteristic  of  a damped  harmonic  oscillator. 

The  detailed  description  of  the  individual  contributions  to  the  real  and  imaginary  parts  of 
the  dielectric  constant  is  given  in  the  Reference  [9].  The  presented  model  is  applicable  in  the 
photon  energy  range  from  0.4eV  to  about  4eV.  The  initial  crystal  parameters  required  for  the 
calculations  include  energy  band  gap  values,  effective  electron  and  hole  masses,  static  dielectric 
constant,  and  the  values  of  the  bowing  parameters  used  for  the  energy  band  gap  calculations  for 
the  ternary  alloys.  The  complex  dielectric  constant  and  refractive  index  of  binary  alloys  (AlAs, 
GaAs,  InAs)  were  first  calculated  and  the  results  were  then  used  in  the  calculation  for  AlGalnAs 
quaternary  alloy.  The  reported  model  also  incorporates  the  lifetime  broadening  effects  of  the  free 
electron-hole  pair  states  through  the  damping  parameters  [12]. 

The  refractive  index  of  a semiconductor,  n(E),  and  the  extinction  coefficient,  k(E),  were 
calculated  in  terms  of  the  complex  dielectric  function  as  follows: 


There  are  a total  of  12  unknown  parameters  required  by  the  model  in  Eqs.(l)  and  (2)  for 
the  real  and  imaginary  part  of  the  dielectric  constant  of  a binary  alloy,  which  include  strength 
and  damping  parameters  for  various  transitions.  These  parameters  are  obtained  by  fitting  through 
minimization  of  the  total  error  function,  F,  over  all  available  experimental  points  for  each 
semiconductor  alloy  of  interest,  which  is  given  below  as: 


i=l 


CALCULATION  RESULTS  AND  DISCUSSION 

The  present  investigation  of  the  optical  properties  of  the  III-V  semiconductor  alloys  is 
centered  on  describing  the  behavior  of  the  refractive  index  for  the  photon  energies  in  the 
transparent  region  as  well  as  for  the  higher  energies.  The  calculated  index  of  refraction  from  our 
extended  ITC  model  and  the  available  experimental  data  for  AlGalnAs  quaternary 
semiconductor  alloy  are  shown  in  Figures  1 and  2 demonstrating  excellent  agreement  between 
the  two.  The  strongest  resonance  peak  of  the  index  of  refraction  occurs  at  the  Ei  transition 
energy.  The  Ei  peaks  in  this  calculation  appear  to  be  generally  sharper  than  the  experimental 
values,  which  can  be  explained  by  the  absence  of  the  lifetime  broadening  effects  in  the  current 
model  for  the  Ei  optical  transition. 

It  is  believed  that  by  adding  the  broadening  constant  as  well  as  including  the  exciton 
effects  for  the  Ei  transition  energy,  better  agreement  for  this  peak  can  be  obtained.  The  increase 
in  the  broadening  parameter  of  AlGalnAs  alloy  is  reported,  which  may  be  attributed  to  the 
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potential  fluctuations  resulting  from  random  atomic  placement  in  the  quaternary  alloys  as 
compared  with  binary  alloys.  Maximum  anion  disorder  occurs  at  y=0.5,  while  that  for  the  cation 
occurs  at  x=0.5  and  y=l  [13].  Therefore,  the  high  atomic  disorder  in  the  triple-cation  sublattice 
of  Alo.3Gao.i6lno.54As  alloy  with  compositions  x,  y,  z close  to  the  middle  point  may  be 
responsible  for  the  slight  deviations  of  the  calculated  refractive  index  from  the  one  determined 
experimentally.  In  this  case  the  bowing  parameters  can  no  longer  be  represented  as  a quadratic 
function  of  the  alloy  composition.  The  atomic  disorder  and  composition  fluctuations  are 
expected  to  broaden  the  optical  spectra. 

The  valley  contributions  to  the  real  part  of  the  dielectric  constant  are  separated  according 
to  the  partition  of  the  Brillouin  zone.  The  L region  with  corresponding  direct  E|  and  indirect  El 
optical  transition  energy  contributes  approximately  65  -75%  to  the  total  value  of  the  dielectric 
constant  while  the  F region  corresponding  to  the  Eo  transition  accounts  to  about  5-10%,  the  X 
region,  which  is  represented  by  indirect  Ex  energy  and  constant  M,  contributes  about  15  to  30% 
of  the  total.  Thus,  the  index  of  refraction  is  essentially  determined  by  the  band  structure  away 
from  the  center  of  the  BZ,  and  modifications  of  the  electronic  structure  at  L and  X points  rather 
than  T CPs,  produce  the  observed  variations  in  the  index  of  refraction. 


Figure  I.  Refractive  index  of  AhGavlnzAs  for  x=0.48,  y=0,  z=0.52. 
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photon  energy,  eV 


Figure  2.  Refractive  index  of  AlxGayInzAs  for  x=0,  y=0.47,  z=0.53. 

The  contribution  to  the  total  dielectric  function  of  3D  discrete  exciton  transitions  close  to 
the  Eo  transition  energy  are  found  to  be  negligible,  about  0.1%,  due  to  the  very  narrow  spectral 
range  of  such  transitions. 

The  results  of  the  AlxGayInzAs  alloy  investigation  were  extrapolated  to  the  different 
compositions  of  interest  to  photonic  devices.  In  order  to  select  the  optimal  materials  for  the 
semiconductor  Distributed  Bragg  Reflectors  (DBRs)  for  application  in  long  wavelength  Vertical 
Cavity  Surface  Emitting  Lasers  (VCSELs),  the  maximum  and  minimum  refractive  indices  were 
calculated  for  response  to  the  incident  photon  energy  of  0.8eV  for  the  alloys  lattice  matched  to 
the  InP  substrate  under  the  direct  band  gap  conditions.  The  Alo.05Gao.42Ino.53As/InP  material 
system  with  index  of  refraction  difference  of  0.46  is  therefore  recommended  for  the  Bragg  mirror 
applications. 

SUMMARY 

An  extended  ITC  model  and  calculations  of  optical  properties  for  AlGalnAs  quaternary 
III-V  semiconductor  alloy  are  presented  and  the  results  are  compared  with  the  experimental  data. 
The  successful  fit  of  the  refractive  indices  in  the  transparent  optical  region  as  well  as  for  the 
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higher  photon  energies  was  attained  by  combining  several  interband  transition  contributions.  The 
largest  contribution  to  the  dielectric  function  is  due  to  the  direct  and  indirect  optical  transitions 
along  <1 1 1>  and  <I00>  directions  in  the  BZ,  which  accounts  for  85-90%  of  the  total 
contributions.  Therefore,  except  for  the  optical  absorption  in  the  vicinity  of  the  T gap,  most  of 
the  optical  properties  of  the  material,  especially  the  index  of  refraction,  are  determined  by  the 
electronic  structure  around  L point,  rather  than  at  the  center  of  the  BZ.  Since  our  model  is  more 
sensitive  to  the  indirect  band  gap  transitions,  the  nature  of  the  indirect  band  gap  contributions  to 
the  total  dielectric  function  of  the  alloy  may  be  better  understood.  The  influence  of  the  discrete 
exciton  states  around  the  Eo  edge  at  room  temperature  was  found  negligible.  We  have  applied  the 
results  of  these  calculations  to  the  AIGalnAs/InP  material  system  for  applications  as 
semiconductor  distributed  Bragg  reflectors. 
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